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CHAPTER  i 

INTRODUCTION 

High  speed  digital  circuits  have  resulted  from  the  miniaturization  of  solid-state  devices 
and  the  implementation  of  those  circuits  into  dense  systems  (VLSI  and  VHS1).  Switching 
devices  presently  used  have  rise  times  and  pulse  widths  in  the  range  of  a  few  nanoseconds. 
Moreover,  devices  with  switching  delays  of  less  than  10  picoseconds  are  being  developed 
from  which  delays  of  less  than  50  picoseconds  will  probably  result  in  large  computers. 
Aside  from  the  physical  characteristics  of  the  devices,  interconnections  are  another  parame¬ 
ter  which  determine  the  speed  of  these  systems.  Today,  interconnections  represent  the  main 
limiting  factor  in  the  reduction  of  speed  because  of  the  inevitable  RC  delay  they  introduce. 
Moreover,  with  increased  density  and  speed,  the  problem  of  cross-coupling  has  become  more 
critical  since  noise  resulting  from  unwanted  signals  must  be  limited.  Finally,  reflections  at 
the  terminals  due  to  mismatch  along  with  device  nonlinearities,  contribute  to  the  complica¬ 
tions  which  have  posed  some  severe  constraints  for  switching  network  design. 

In  this  study,  the  problems  associated  w’ith  the  propagation  of  fast  pulses  in  microstrip 
interconnections  will  be  analyzed.  The  purpose  of  this  analysis  is  to  investigate  each  effect 
such  as  RC  delay,  reflections,  cross-coupling,  nonlinearities,and  to  propose  a  model  for  the 
combined  effect.  Such  a  model  is  essential  since  it  provides  the  engineer  with  the  necessary 
elements  from  which  design  guidelines  can  be  implemented  while  giving  a  better  under¬ 
standing  of  noise  phenomenon. 

In  order  to  develop  the  model,  it  will  be  necessary  to  determine  the  physical  characteris¬ 
tics  of  the  interconnections  as  well  as  those  of  the  propagation  medium.  Unfortunately, 
analytical  expressions  predicting  those  physical  characteristics  are  not  always  realistic  and 
do  not  account  for  processing  variations,  nor  can  they  be  applied  to  arbitrary  geometries.  To 


circumvent  this  obstacle,  several  measurement  techniques  will  have  to  be  developed  in  order 
to  obtain  the  important  physical  parameters. 

An  overview  of  the  presently  available  microstrip  technology  along  with  the  theoreti¬ 
cal  foundations  is  first  explained.  Some  techniques  of  measurements  in  the  time  domain  and 
in  the  frequency  domain  are  then  discussed  for  both  single  and  coupled  lines.  Experimental 
results  are  finally  analyzed  and  extensions  to  multiple  lines  are  attempted. 

1.1  Characteristics  of  Microtrip  Transmission  Lines 

Numerous  investigators  have  studied  the  properties  of  microstrip  transmission  lines  ,  pro¬ 
posed  models,  design  equations,  and  various  approaches  for  the  analysis.  Figure  1.1  shows  a 
cross  section  of  a  microstrip  line  which  consists  of  a  strip  conductor  of  width  W  resting  on 
a  substrate  of  thickness  h  and  dielectric  constant  e, ,  and  a  ground  plane  which  serves  as  the 
return  path.  The  conductor  can  be  either  copper,  gold,  or  an  alloy  of  tin  and  nickel. 

Since  the  fields  between  the  strip  and  the  ground  plane  are  not  entirely  within  the  sub¬ 
strate,  the  propagating  mode  along  the  strip  is  not  purely  tranverse  electromagnectic  (TEM) 
but  quasi-TEM  with  a  phase  velocity  v0  given  by 


where  c  is  the  speed  of  light  in  vacuum  and  e,, ,  the  effective  dielectric  constant  of  the  sub¬ 
strate.  The  wavelength  X„  in  the  line  is  : 


(1.2) 


where  /  is  the  frequency  of  the  propagating  signal.  Moreover,  if  the  line  is  assumed  to  be 
lossless  its  characteristic  impedance  (or  wave  impedance)  is  expressed  by 


4 


(ohms) 


W/h 


Figure  1.2  :  Characteristic  impedance  /„  versus  width  to  height  ;  IV  'h.  Solid  line  :  theoreti¬ 
cal  expressions  from  ilammerstud  with  €  =4.5  ;  points:  experimental  measure¬ 
ment  using  time  domain  reflectometrv.  l  or  the  experimental  points,  the  lines 
hace  a  dielectric  constant  4$e  ^5. 
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For  W  /h  >  2, 


^ -  =  —  [B  ~  1  -  ln(2£  -  1)  +  — - — -(ln(5  -  1)  +  0.39  -  ^-)] 
h  v  2e,  €, 


(1.6) 


where 


A 


€,  -1 

€,+ 1 


(0.23  -  -2^1) 

€, 


(1.7) 


and 


_  377tt 


(1.8) 


These  relations  assume  an  infinitely  thin  conductor  strip  ;  in  practice,  the  thickness  must  be 
considered,  which  is  done  by  defining  an  effective  width  of  the  strip,  VV,  for  use  with  Equa¬ 
tions  (1.5)  -(1.8)  instead  of  W . 


For  IV  /h  -±-, 
2v 


Wc  =  +  -^-(1  +  2-) 

h  Trh  t 


(1.9) 


For  W  /h  < 

27 T 


+  -Ml  + 

nh  t 


(1.10) 


It  has  been  found  that  for  t /h  <  0.005,2  <  €,  *5  1  .and  0.1  /h  ^  0.5  ,  neglecting 


the  thickness  of  the  strip  introduced  very  little  error.  Since  these  values  represent  most 
practical  cases,  correction  for  finite  thickness  is  often  neglected. 

1.3  Frequency  Dependence  and  Dispersion 

The  relations  derived  above  were  based  on  the  quasi-TEM  model.  At  higher  frequencies 
this  model  is  not  accurate  for  microstrip  lines  because  of  the  propagation  of  hybrid  modes. 
As  the  frequency  increases,  the  phase  velocity  decreases  which  reveals  an  increase  of  the 
effective  dielectric  constant ;  an  increase  of  the  characteristic  impedance  is  also  observed. 

The  frequency  /„  below  which  dispersion  effects  can  be  ignored  is  given  by  fll] 


f0(Ghz )  =  0.3 

V  Mcm)VT-l 

and  the  effective  dielectric  constant  €,,(/  ),  is  given  by  [l 2] 
e,  — 

€,  =  €, - 

1  +C( 

where 


(1.11) 


(1.12a) 


(1.12b) 


G  =  0.6  +  0.009Z, 


(1.12c) 


The  frequency  dependence  of  the  characteristic  impedance  Z„,  was  reported  by  Owens  [13] 
y(f)  — _ yiih _ 

’ 1  111 


where  VVt. (/  )  is  the  effective  width  given  by 


1.4  Time  and  Frequency  Domain  Characterization 

The  steady-state,  single  frequency  solutions  for  the  voltage  and  current  on  a  transmission 
line  derived  from  Maxwell’s  equations  [14]  have  the  form 


v(tjc)  =  [A  e  '°  +  b  e  '°]e’ 


/(/,*)= -J-U  e  '’-Be  ,"]eJ0' 


(1.15) 


(1.16) 


where  A  and  B  are  constants  which  depend  on  the  terminations.  This  can  be  rewritten  in  a 
more  general  form  which  accounts  for  the  transients  in  the  line  : 


/  .  \  i r  /  .  \  .  *  >  /  .  \ 


3  :  Time  Domain  Retlectometrv.  Applications  to  transmission  line 

characterization^)  experimental  setup  ,  (b)  voltage  reading.  The  time  axis  can 
also  be  used  as  position  axis  which  allows  to  determine  and  locate  discontinuities 
along  the  test  line.  Rise  time  of  applied  step  ==50  ps  . 


ramp  generator  V(  of  internal  resistance  Zg  (usually  50  ft).  A  coaxial  cable  of  characteris¬ 
tic  impedance  Z,.  and  length  l  is  connected  to  the  line  to  be  tested  (characteristic  impedance 
ZQ  and  length  1).  When  the  switch  is  closed  at  t  =t0,  a  forward  moving  wave  C  +  is  gen¬ 
erated  at  point  a .  For  voltage  and  current  continuity,  it  must  satisfy 


V,  -c+  =  C+ 
zt  Z, 

from  which  we  obtain 


(1.19) 


(1.20) 


Upon  reaching  the  test  line  at  point  b,  a  backward  wave  C_  as  well  as  a  forward  wave  T  + 
appear.  They  must  satisfy 


C+  +  C_  =  7\ 


(1.21a) 


(1.21b) 


which  yields 


C_ 


=  ( 


20  -Z, 

Z0  +  zt 


)c+  =  r,c  + 


(1.22a) 


7\ 


2  Z„C  + 
7.,,  +  Z, 


(1.22b) 


where  r„  is  the  reflection  coefficient  from  the  test  line  at  b.  When  l\  reaches  the  end  of  the 


test  line  at  point  c,  a  reflected  wave  T_  is  produced  such  that 

r_  =  (^  yr+  =  rLT+  (1.23) 

*-L  +  Zo 

where  Fz  is  the  reflection  coefficient  at  the  load.  Upon  reaching  b,T _  will  generate  two 
new  waves  :  C_  on  the  coaxial  line  and  7”+on  the  test  line.  These  successive  reflections  will 
keep  taking  place  indefinitely.  The  expression  for  the  voltage  at  any  point  along  the  line 
w'ill  have  the  form  of  a  geometric  series  of  T,  rL  which  converges  to  the  final  value  (for 
very  large  t) 


(1.24) 


Figure  1.3b  shows  a  graph  of  the  voltage  at  point  a,  Va(t )  as  a  function  of  time.  Since  the 
phase  velocity  in  both  coaxial  and  test  lines  are  assumed  constant,  there  exist  a  linear  rela¬ 
tionship  between  time  and  position.  The  first  step  in  the  voltage  is  the  magnitude  of  the  for¬ 
ward  w'ave  traveling  the  coaxial  cable  C+.  Its  duration  is  2 l/vc  where  v(  is  the  velocity  in 
the  coaxial  line.  The  second  step  is  T +.  If  the  magnitudes  are  normalized  to  C+,  the 
difference  between  C+  and  T  +  is  the  reflection  coefficient  Tn  due  to  the  test  line 


4,  -  Z„ 
+  Z, 


(1.25) 


When  f„  is  read  from  the  Time  Domain  Reflectometer.  it  allows  direct  calculation  of  Z„  : 
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Whereas  the  phase  velocity  va  can  be  obtained  from  the  time  delay  t„  : 


(1.27) 


Capacitance  and  inductance  per  unit  length  can  be  extracted  using  (1.1  Ml. 3). 


(1.28a) 


L  m*L 

'  vn 


(1.28b) 


1.4.2  Frequency  Domain  -  Scattering  Parameter  Model 

The  information  obtained  from  the  TDR,  although  complete,  does  not  provide  us  with 
any  insight  of  the  frequency  dependence  of  the  transmission  line  parameters.  Because  of  ins¬ 
tability  concerns  at  large  f requencies  (RF  and  microwave)  and  because  of  their  manoeuvra¬ 
bility  at  those  frequencies.  Scattering  parameters  (S  parameters)  can  be  used  for  a  frequency 
domain  analysis.  A  good  introduction  to  S  parameters  is  given  in  References  [16]  and  [17], 
We  briefly  summarize  the  essentials  of  these  parameters  and  their  meanings. 

Figure  1.4  shows  a  two-port  network  inserted  on  a  line  of  characteristic  impedance  Z„. 
Both  ends  of  the  line  can  be  regarded  as  connected  to  generators  of  internal  impedance  Z. . 
E, ,  and  E,2  are  the  incident  portions  of  the  voltages  in  ports  1  and  2,  respectively.  Simi¬ 
larly,  E,  ,  and  E,  2  are  the  reflected  portions  of  the  same  voltages.  We  can  then  define  four 
traveling  waves  on  the  line  : 


vV\V- .  -Y 


V-*  ’/-V  .  ,N  .  ' 


PORT  I 


PORT  2 


bi  =  ^hqi  +  22a2 

b2S  ^l0!  +  S22°2 


Figure  1.4  :  Scattering  matrix.  Definition  in  terms  of  the  power  waves. 


zQ-z 

Z0+Z 


Figure  1.5  .-Experimental  setup  for  determining  the  characteristic  of  a  microstnp  line  as  a 
function  ol  frequency.  Equations  are  shown  that  relate  the  velocity  and  the 
impedance  to  the  measured  reflection  and  transmission  parameters. 


The  square  of  the  magnitude  of  these  new  variables  has  the  dimension  of  power  ;  I  a,  1 1  is 
the  incident  power  to  port  1  ;  I  b  2 1 2  is  the  reflected  power  f rom  port  2  of  the  network.  The 
four  S  parameters  characterizing  the  two-port  relate  the  traveling  waves  as  follows  : 


b  |  =  S  ||U  i  +  S  2  (1.30a) 

b2  =  S  zia  i  +  S2**2  ( 1.30b) 

From  these  relations,  the  four  S  parameters  can  be  defined  alternatively  as 


(1.30c) 


by  b , 

Si2=—  S  ->2  —  ~  (1.30d) 

0-2  a  i=0  a  2  a  i=0 

If  both  ends  of  the  line  are  connected  to  generators  with  internal  impedance  Ze ,  then  a  ,  and 
a  2  are  simply  the  incident  waves  from  the  generators.  Since  the  incident  power  from  either 
generator  can  be  deliberately  set  to  zero.  Equations  (1.30c)  and  (1.30d)  provide  the  main 
scheme  for  determining  the  S  parameters  of  a  network.  It  is  of  interest  to  note  that  S ,,  and 
S  22  are  reflection  coefficients  and  relate  traveling  waves  of  one  given  port  whereas  S  l2  and 
S  21  are  the  transmission  coefficients  of  the  two-port. 


1.4.3  Scattering  Matrix  for  Transmission  Line 

When  the  network  to  be  characterized  is  a  transmission  line,  it  can  be  modeled  as  a  two- 
port  (Fig.  1.5).  Consider  an  unknown  transmission  line  of  characteristic  impedance  Z0  and 
length  l  embedded  between  two  coaxial  lines  of  impedance  Zg  terminated  with  their 
characteristic  impedance.  On  the  left  is  an  RF  source  of  frequency  /  =<u/27 r.  Setting  the  ori¬ 
gin  at  the  first  intersection  ,  the  voltage  and  current  equations  for  the  system  are 


where 
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u / 


0  =  e  • 


and 


(1.34a) 


T  =  (1.34b) 

Z0  +  Zg 

r  is  the  reflection  coefficient  at  the  ends  of  the  unknown  line.  Adding  (1.33)  to  (1.34),  we 
get 

=  (U5) 
Subtracting  (1.34)  from  (1.33)  gives 


S2i-Su  =  -f—^  (1.36) 

Combining  (1.35)  and  (1.36)  yields 


(1.37a) 

(1.37b) 


Since  S  u  and  S  i,  can  be  measured  with  a  network  analyzer.  Equation  (1.37)  allows  a  direct 
determination  of  T  from  which  Z„  can  be  obtained.  The  propagation  velocity  and  the 
effective  dielectric  constant  can  also  be  obtained.  Figure  1.6  shows  experimental  plots  of 
characteristic  impedance  versus  frequency  for  various  geometries  of  lines.  The  measure¬ 
ments  were  performed  on  an  HP  8505  automated  network  analyzer.  At  higher  frequencies. 


Z0  (ohm) 


,X  W/h  =  0  61 


•o  W/h  =  0.92 


-a  W/h  =  128 


FREQUENCY  (GHz) 


f  igure  1.6  :  Experimental  plot  of  the  characteristic  impedance  Z„  as  a  function  of  frequencv 
for  different  geometries.  The  dashed  tines  are  the  best  fit  curves.  Measurements 
were  performed  on  an  HP  8505  RF  network  analyzer  and  Equation  (1.37) 
allowed  determination  of  Z,  . 


the  measurements  , become  more  difficult  and  external  effects  such  as  junction  discontinuities, 
capacitance,  and  inductance  become  more  significant. 
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CHAPTER  2 


COUPLED  TRANSMISSION  LINE  CHARACTERISTICS 


The  theory  of  coupled  transmission  lines  arose  with  the  early  applications  to  multicon¬ 
ductor  systems.  Its  development  continued  with  the  advent  of  directional  electromagnetic 
couplers  [18].  Several  methods  using  the  quasi-TEM  approximation  have  been  proposed  to 
obtain  the  parameters  of  coupled  microstrip  lines.  A  comparison  of  these  methods  is  given  in 
[19].  Krage  and  Haddad  [20],  [21]  have  determined  the  inductive  and  capacitive  coupling 
coefficients  as  well  as  the  directivity  for  various  geometries  of  coupled  lines.  Bryant  and 
Weiss  [22]  have  established  the  relationship  between  the  electrical  and  physical  parameters 
of  coupled  line  pairs  using  the  Green’s  function  approach.  Their  .V1STRIP  computer  program 
has  been  validated  after  comparison  with  various  authors.  Garg  and  Bahl  [23]  derived  semi- 
empirical  equations  for  the  even  and  odd  mode  parameters  Hammerstad  and  Jensen  [24]  suc¬ 
ceeded  in  implementing  a  model  with  errors  iess  than  those  caused  by  physical  tolerance. 
Recently,  kirshing  and  Jensen  [25]  reported  frequency  dependent  expressions  with  unpre¬ 
cedented  accuracy.  Most  of  the  solutions  provided  by  these  numerical  techniques  suffer  from 
the  lack  of  experimental  data  which  would  validate  their  application  to  microstrip  coupler 
design.  The  main  difficulty  arises  with  the  presence  of  different  modes  of  propagation  and 
the  parameters  associated  with  these  modes  which  impose  more  complex  measurement  tech¬ 
nique  requirements.  The  goal  of  this  chapter  consists  of  implementing  the  electrical  model 
for  coupled  lines  and  determining  experimentally  the  coupling  parameters.  Correlations 
between  physical  dimensions  and  electrical  parameters  can  then  be  established  empirically. 
Altough  this  method  limits  the  flexibility  of  the  microstrip  designer,  it  provides  more  accu¬ 
racy  for  parameter  determination  while  giving  a  better  insight  of  the  coupling  phenomena. 
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(b)  CIRCUIT  MODEL 


Figure  2.1  .  Coupled  microstrip  lines  :  (a)  geometry,  (b)  distributed  circuit  equivalent.  It  is 
assumed  that  the  two  lines  are  identical.  /..  and  C4  are  the  self-inductance  and 
capacitance  of  each  line.  Ln,  and  C,„  are  the  mutual  parameters  of  the  pair 
These  parameters  can  be  measured  and  related  to  the  matrix  elements  usin^ 
Equations  12.5  W  2.8).  s 


Comparison  with  Equations  (2.1M2.2)  with  (2.3M2.4)  gives  the  equivalences. 


C  n  —  Cs  +  Cm 

(2.5) 

c  12  =  —C,n 

(2.6) 

II 

(2.7) 

t- 

II 

t- 

(2.8) 

For  the  symmetrical  case  one  solution  can 

be  found  by  adding  (2.1a)  to  (2.1b)  and  (2.2a)  to 

(2.2b). 

ex  ■  si 

(2.9a) 

=<c„+Ci;)31'- 

dx  at 

(2.9b) 

where 

li 

+ 

< 

(2.10a) 

/.,  —  ~(/  |  +  /  ■*) 


(2.10b) 
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are  the  common  mode  voltage  and  current.  Equations  (2.9a)  and  (2.9b)  have  the  same  form 
as  the  well-known  single  line  Telegraph  Equations  ;  therefore  the  solutions  in  the  fre¬ 
quency  domain  assuming  a  time  harmonic  dependence  are 


<4>.r  ^  .  <i>  r 

ve(x)  =  ac  e  ‘‘  +  Bce  '' 


(2.11a) 


(2.11b) 


where  Ac  and  Bc  are  constant  coefficients  associated  with  the  forward  and  backward  trav¬ 
eling  waves,  respectively.  Z0  and  v(,  are  the  even  mode  impedance  and  propagation  velocity, 
respectively  : 


(2.12) 


KL  i,+L  12XC  „+C  ,2)  JU~*LJC~ 


(2.13) 


and  cu=27 r/  is  the  angular  frequency  of  the  propagating  signal.  Subtracting  (2.1b)  from 
(2.1a)  gives 


dx 


=  an-Lu)^- 


(2.14a) 


Mj. 

dx 


(2.14b) 


where 


are  the  differential  mode  voltage  and  current.  As  before  this  system  can  be  solved  to  give 


vAx)  =  Ad  e  '*  +  Bd  e  J"d 


(2.16a] 


A  ~J —  U  — 


(2.16b) 


Ad  and  Bd  are  as  before  constant  coefficients  associated  with  the  forward  and  backward 
waves,  respectively,  and  Zd  and  vd  are  the  odd  mode  impedance  and  propagation  velocity 
defined  as 


y/(L  n-L  l2XC  n-C  l2) 

—Lm  XC,  +2Cn, )  U,18; 

Using  Equations  (2.10)  and  (2.15) 

we  obtain  the  general  solution  for  the  line  currents  and 

voltages  : 

V,(x)*  V,U)  +  Vd(x ) 

(2.19a) 

/  |(x  )  =  /p(x  )  +  Id(x  ) 

(2.19b) 
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Vlx)  =  Vc(x)-  Vd(x) 

(2.20a) 

H, 

1 

>< 

II 

--“v 

H 

(2.20b) 

which  gives 

v  ,(* )  =  ac  e  +  B'  e  +  Ad  e  "d  +  Bd  e 


+j— 


(2.21a) 


(2.21b) 


V  U  )  =  Ac 


+j- 


+  B,  e 


-  a,  e 


tux 

/ 

vd 


—  bh 


(2.22a) 


(2.22b) 


Equations  (2.21)  and  (2.22)  give  the  general  form  for  the  line  currents  and  voltages  for  a 
time  harmonic  excitation.  The  constants  Af,  Be,  Ad ,  and  Bd  are  determined  by  the  boun¬ 
dary  conditions  at  the  four  ends  of  the  lines.  Determination  of  these  constants  depends  on 
the  nature  of  the  problem.  For  the  case  where  the  excitation  is  periodic,  a  Fourier  series  form 
must  be  assumed  as  the  general  solution  for  which  the  constants  become  the  Fourier  series 
coefficients.  Finally  for  nonperiodic  excitation,  a  Fourier  transform  approach  must  be  used 
for  the  general  solution.  This  will  be  illustrated  in  the  next  chapter.  In  all  cases,  the  time- 
dependent  solutions  have  the  form  : 


v,(x)  =  aPu  -  -£-)+  bjt  +  — )  +  adu  -  — )  +  bdu  +  — ) 

V,,  Vj  vd 


(2.23a) 


?S 


*?T  j*j"j*r  v*W  fj  u"J'*m  u  m  •»•:■■-'  *  t*  m  ^tr*  it  *  t*  TT^rf  r* 


‘J*\ 


R 


>  V. 

k' 


I  t(x)  =  y-lajt 
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—)-b,Ai  +-*-))  + J-[dAt 

Vc  Vc  Zj 


-JL)-bAt  +  JLj] 


vix)=acu  -  —)+bcu  +  —)-adu  -  -*-)-bd(t  +  — ) 


(2.23b) 


(2.24a) 


/ A )  =  a  6Pa  +  -2-)] - ^-[a,u  6rf(*  +  — )]  (2.24b) 

Zc  vf  v„  Zd  vd  17  vd 

CLc  and  CLd  are  the  forward  functions  for  the  even  and  odd  modes,  respectively,  6C  and  bd 
are  the  backward  functions  for  the  same  modes. 


2.2  Application  and  Coupling  Parameter  Measurement 

The  major  advantage  of  Equations  (2.23)  and  (2.24)  is  that  the  solutions  are  put  in  a 
form  where  the  reflected  waves  are  separated  in  time  from  the  forward  waves.  This  can  be 
applied  w  ith  TDR  measurements  to  obtain  the  coupling  parameters  of  a  microstrip  pair  since 
only  the  forward  voltage  waves  need  to  be  considered.  In  Figure  2.2a,  two  identical  coupled 
lines  of  length  l  are  connected  at  x  =  0.  A  coaxial  cable  of  characteristic  impedance  Zg  con¬ 
nects  the  TDR  to  the  pair.  A  forward  traveling  voltage  step  of  magnitude  V  f  upon  reach¬ 
ing  the  junction  generates  a  backward  step,  Vb.  From  Kirchhoff's  law  we  must  have 


Vtd,  =  <20(f,0)  +  ad(t,0)  =  V  j  +  Vs 


*  xd t  — 

a,  +  a. 

+ 

a,  <•  -0)  a, <•  •0) 

r  •: 

zr  “  z„ 

c  7  «  7 

*~d 

i*ts 

T-' 

■  V,J,  = 

CLc  (t  ,0)  —  ad(t, 0) 

V,  -V* 


(2.25a) 


(2.25b) 


(2.26) 


j  -  jiLiX  a— »  ■  «■*- 


This  implies  that 
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Ze  *  2 


ve  *  time  delay 
(NO  ODD  MODE) 


(a)  EVEN  MODE  IMPEDANCE  DETERMINATION 


Ze  Zd 
Ze+Zd 


(b)  "parallel"  CONNECTION 


(i±£z.)z 


Figure  2.2  :  Experimental  determination  of  the  even  and  odd  mode  characteristics  of  a  pair 
of  coupled  transmission  lines  :  (a)  even  mode  measurement  since  no  odd  mode  is 
present  both  velocity  and  impedance  for  this  mode  can  be  measured  :  <fr> parallel 
connection  since  it  is  difficult  to  excite  the  odd  mode,  this  configuration  is  used  to 
determine  the  parallel  combination  of  Z„m  and  Zrjdd. 


aAt. o)  =  o 


U.//J 


and 


V,«, 

I  idir 


2, 

2 


(2.28) 


Therefore,  the  reflection  coefficient  measured  at  x  =0,  pc  =Vh /V  f  can  be  related  to  the  even 
mode  impedance  of  the  pair  : 


Z,  =  2(^-)Ze  (2.29) 

1  ~Pe 

Since  no  odd  mode  exists,  the  measured  time  delay  tc,  can  be  related  to  the  even  mode  velo¬ 
city  : 


(2.30) 


Another  configuration  is  that  of  Figure  2.2b  in  which  line  2  is  connected  to  ground  at  x  *  0 
and  line  1  is  connected  to  the  coaxial  line.  The  equations  are 


v,d>  =  a,  o.o)  +  ad(t,o)  =  vf  +  v,. 


(2.31a) 
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=  Z'Zj 
1  Id.  ZK.  Zj 


(2.32b) 


The  associated  reflection  coefficient  p;,  =  Vh/V }  can  be  measured  directly  and  using  the 
above  relations  we  obtain 


ZeZj  _ 
Z..+2,  * 


1+gg 

*-Pp 


(2.33) 


Since  both  even  and  odd  modes  are  propagating  ,  a  direct  relation  between  the  time  delay 
and  either  velocity  cannot  be  established.  Nevertheless,  Equations  (2.29),  (2.30),  and  (2.33) 
provide  the  necessary  relationship  for  extracting  Zc,  Zd,  v,,  and,  if  the  single  line  parame¬ 
ters  are  known.  Equations  (2.12),  (2.13),  (2.17),  and  (2.18)  can  be  used  to  determine  the 
mutual  capacitance  and  inductance  as  well  as  the  odd  mode  propagation  velocity  of  the  pair 
(Figure  2.3).  From  these  equations  we  see  that 


Zd  <  Z„  <  Zr  (2.34) 

which  suggests  that  for  weak  coupling 


Zj  =  Z„  ==  Z„  (2.35) 

Table  1  shows  experimental  values  of  self  and  coupling  parameters  obtained  for  several 
microstrip  pairs.  These  values  were  extracted  by  using  the  methods  outlined  above.  As 
expected  mutual  parameters  depend  strongly  on  the  spacing  between  the  lines  ;  it  is  also 
observed  that  the  mutual  inductance  is  stronger  even  at  wider  spacings  between  the  edges  of 


the  coupled  lines. 


TABLE  1  :  EXPERIMENTAL  DATA  OBTAINED  FOR  VARIOUS  GEOMETRIES  OF 
TRANSMISSION  LINES  WITH  GLASS-EPOXY  DIELECTRIC.  GEOMETRICAL 
DIMENSIONS  ARE  AS  DEFINED  IN  FIGURE  2.1,  THEIR  UNIT  IS  IN  MILS. 
IMPEDANCES  ARE  IN  OHMS  AND  VELOCITIES  ARE  IN  m  /ns.  INDUC¬ 
TANCES  AND  CAPACITANCES  ARE  IN  nH  lm  AND  pF  lm ,  RESPECTIVELY. 
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Since  different  charge  distributions  arise  with  the  two  modes,  some  authors  have  found  it 
convenient  to  define  an  even  mode  capacitance  and  C,.,  and  an  odd  mode  capaciance  CtMj  and 
base  the  complete  analysis  in  terms  of  these  parameters  [8],  [9].  Figure  2.4  shows  a  decompo¬ 
sition  of  the  total  capacitance  of  a  microstrip  pair  in  terms  of  even  and  odd  mode  capaci¬ 
tances.  The  existing  expressions  for  C„  and  C ^  are  semi-empirical,  however,  using  our 
model,  it  can  be  shown  that  C...  -  Cs  and  Crjd  -  C,  +  2 Cm. 

The  major  advantage  of  the  techniques  of  measurement  introduced  above  is  that  they  are 
very  accurate,  since  no  approximation  was  made  in  deriving  the  model  and  the  relations 
between  the  coupling  and  propagation  parameters.  As  a  consequence,  reliable  empirical  data 
can  be  established  from  these  measurements  and  used  as  design  guidelines  for  microstnp 
couplers. 


(a)  EVEN  MODE 


(b)  ODD  MODE 

Figure  2.4  :  Decomposition  of  the  capacitances  into  odd  and  even  mode  components.  Cr 
represents  the  capacitance  associated  with  the  TFM  propagation  mode  .  C.  is  the 
fringing  capacitance.  C?0  and  C..d  are  the  mutual  capacitances  in  air  and  in  the 
dielectric,  respectively.  The  normalized  charge  distribution  associated  with  each 
mode  is  shown  on  the  right  with  <r„  being  the  charge  at  the  center  of  each  strip. 
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CHAPTER  3 

TRANSIENTS  IN  COUPLED  MICROSTRIP  LINES 

In  the  preceding  chapter,  the  response  of  coupled  lines  was  analyzed  for  harmonic  (single 
frequency)  excitation.  We  now  wish  to  extend  the  analysis  to  account  for  the  transients 
involved  in  the  case  of  an  arbitrary  excitation.  This  is  particularly  important  for  digital 
applications  where  the  signal  applied  to  microstrip  lines  is  a  pulse  or  an  impulse.  Since  the 
rise  and  fall  times  of  these  excitations  are  becoming  shorter,  the  coupling  between  adjacent 
lines  becomes  more  significant  and  has  more  serious  effects.  This  is  the  so-called  "crosstalk" 
noise.  Various  authors  have  attempted  to  describe  this  coupling  phenomenon  and  relate  it  to 
the  mutual  parameters  of  the  pair.  Cotte  [26]  developed  a  first  theory  on  the  propagation  of 
pulses  in  a  coupled  pair  of  conductors.  Catt  [27]  analyzed  the  same  effects  for  various 
transmission  line  configurations.  Jarvis  [28]  studied  the  waveform  distortion  caused  by  this 
coupling  as  well  as  the  effects  of  the  terminations.  Several  numerical  techniques  have  also 
been  applied  for  simulation  of  these  coupling  effects.  In  general,  the  cross-coupling  between 
two  lines  is  a  function  of  the  terminations  which  make  up  the  boundary  conditions  for  the 
general  solution  previously  derived.  This  chapter  will  examine  the  solution  to  this  problem 
in  the  time  domain. 

3.1  Passive  Terminations 

Many  practical  situations  in  digital  applications  involve  pairs  of  micrcstrip  lines  for 
which  the  the  behavior  of  the  terminations  does  not  vary  with  the  magnitude  of  the  sig¬ 
nals.  In  this  case  a  real  linear  impedance  can  be  used  to  model  the  terminations  at  the  source 
and  at  the  load.  Figure  3.1  shows  a  microstrip  pair  having  resistances  Zs ,  and  Z,:  at  the 
sending  end  and  ZL  ,  and  ZLl  at  the  far  end.  The  voltage  sources  have  magnitudes  Vs ,  and 
V,  2-  The  time  harmonic  solution  for  the  voltages  and  currents  was  found  to  be 
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V  Ax)  =  A,  e  ''  +  B ,  <?  ''  +  Av  £  rf  +  BJ  e  'J 


cji  tu.i  u>i  uii 

4  -  D  +  7 A  “7 D  +  7 

/.(x)-^e  ^ 


v/x )  =  Af  e  +  bc  e  -  Ad  e  d  -  Bd  e 


lix>.  ± 

zc  z,  Z^  ztf 


If  we  apply  the  boundary  conditions  at  x  =0  and  at  x  =Z  ,  namely 


vsl  =  V/0)  +  Zsl/,(0) 

vs2  =  v/o)  +  zs2i  A  o) 

0=  V,(Z)-Z£  ,/,(!) 

0  =  V  U)-ZL2I  ll ) 
and  by  setting 


0.,  =  e  a„  =e 


(3.1a) 

(3.1b) 

(3.2a) 

(3.2b) 

(3.3a) 

(3.3b) 

(3.3c) 

(3.3d) 

(3.4a) 

(3.4b) 

<3.4c) 
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ed  =e  rj  ad  =  e 


(3.4d) 


we  arrive  at 


(l+/n)  (l  —  m)  (1+s)  (1  —  s) 

(1  +  v)  (1  -  v)  -(1  +  y)  (y  -  1) 

Or(  1—q)  ac(l  +  q)  Od(l  —  z)  arf(l  +  z) 

0(,(l  —  d)  aff(l  +  d )  9d(p  —  l)  —a^(l  +  p) 


which  is  the  generalized  matrix  equation  for  the  coefficients  Ac,  Be,  Ad,  and  Bd. 

In  many  practical  cases,  the  terminations  on  both  lines  are  identical,  or  Z,  ,=ZS  2=Z,  and 
ZL  \=Zi2~ZL.  This  reduces  Equation  (3.5)  to  a  much  simpler  form.  Then  Cramer’s  Deter¬ 
minant  Rule  can  be  applied  to  give 


(Vs  [  +  Vs2)Z, 


2(zp  +  zs )  i-r„r  Leo* 


(3.6a) 


-<v,l  +  vt2')Zr9r2rL 
2 (Z,  +Z, ) 


i  -r„r,.0‘ 


(3.6b) 


(Vsl-Vs2)ZJ  1 

2{Zd  +  Z, )  \-TidTLd9} 


(3.7a) 


-(vsl- v;2)z,ejr^ 

Zd  +  Zj 


1  -r,„r  Ldej 


(3.7b) 
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r  _  h,  ~Z- 

Lt  ~  ZL  +  Zc 
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Zz  ~Z, 
2/  +  Z,, 


e=e 


od  =e 


(3.8b) 


(3.8c) 


(3.8d) 


Equations  (3.6)  and  (3.7)  give  the  coefficients  for  the  time  harmonic  periodic  excitation 
where  Vs ,  and  Vs2  are  assumed  to  have  the  same  frequency.  From  these  equations  it  is  seen 
that  even  or  odd  mode  can  be  matched.  Particularly,  when  Z,  =Z,  and  ZL  =  Zd,  we  have 


Ar  = 


= 


Aj  — 


(V,,  +  Vs,) 


-<^«  +  vs2wrld 


(vsl-vs2)Zd 


2  (Z*  +Ze) 


Bj  =  0 


(3.9a) 


(3.9b) 


(3.9c) 


(3-9d) 


In  the  case  where  the  excitation  is  not  periodic,  contributions  from  all  frequencies  must  be 
included  in  the  general  solution.  The  coefficients  ,  thus  become  continuous  functions  of  fre¬ 
quency.  In  the  time  domain,  solutions  for  the  voltages  and  currents  are  obtained  by 
integrating  over  all  frequencies.  We  therefore  have 


/+/OJ»r - —  ■  »  -rjufi  - - 

Ae(io)e  >v  da,  +  J  Be(u»)e  l>  da, 


(3.10a) 
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/  Ad(aj) 


+  Jhit  -  —  >  .  +JU*t  +  —I 

6  d  dii>  +  J  Bd(u>)  6  d  d  a> 


1  r  ~  — 1  ,  r  +/«*'  +  — ' 1 

li(tjc)=  -j-J  Ae(u>)e  ’  doi  -~J  Be(io)e  w*  do 


(3.10b) 


i  r  +y“lf  — 1  i  r 

+  —  J  V"  -  —  J  Bd( gj)i 


+ ;<JI  +  - • 


+oo  x  +oo 

.»  +jult - )  .  +  2utr  + — I 

'£tjc)=  J  Ae(a>)e  v'  dui  +  J 


(3.11a) 


+°°  »  +°°  , 
- 1  -  +/<Jr  +  —i 

A^(w)e  *d  dw  -  J  Bd(u>)e  "d  a 


i  r  — 1  .  . 

Af(a>)e  v'  da>  -j~J  Bc(a>)e 


+  yaj»r  + 


(3.11b) 


-j-jAd(a>)e  adw+-±-J  Bd(w)e 


+  yw»r  +  — 1 


This  can  be  rewritten  in  a  simpler  form  as 


V1  =  SJ+S2  +  S,1  +  S,; 


(3.12a) 


S  |  S  2  51 1  S'  2 

zT  ”  zT  +  z7  ~  z7 


(3.12b) 


V,  =  5,  +52-S,1  -  S'; 


(3.13a) 


.  5,  S,  S',  S'2 

2  z,  z„  zd  +  zd 


(3.13b) 
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(3.24b) 


A,  («)  -  2(2^  +  2j )  E  ri„  e 


(3.25a) 


(3.25b) 


Using  Equations  (3.14)  through  (3.17)  and  making  use  of  the  time  delay  theorem  for  Fourier 
transforms,  we  can  then  invert  for  the  associated  modal  voltage  waves  in  the  time  domain. 


5 ,0 , jc )  =  -:(7Zf,7  <  £[vs  ,0  -  £&£)  +  vs it  -  )]  r*  rff 

ZV4eT4j;  *  =o  vc 


(3.26) 


Sitjc)m  -  §IV„(t  +  -■  -■a+l)Z)  +  Vj2(t  +  r*ri+>  (3.27) 

^  =o  Ve  Ve 


S,,(r.*)  = 


2(Zd+Zs) 


£[VS  ,(t  -  -  vs/f  -  ^— -)]  r*,ri 


(3.28) 


?itjc)=  T(~y-<  S[V„(t  +  -  2(k+i)l)-V,it  +  r^rf/1  (3.29) 


and  using  relations  (3.12)  and  (3.13),  the  voltages  and  currents  on  both  lines  can  be  obtained. 

These  equations  show  that  the  resulting  signals  on  the  lines  can  be  expressed  as  an 
infinite  sum  of  delayed  and  attenuated  replicas  of  the  the  original  applied  voltages.  They 
apply  for  arbitrary  waveforms  and  allow  to  determine  the  voltage  and  current  magnitudes 
at  any  time  and  any  position  along  the  lines.  Convergence  of  the  series  depends  on  the 
reflection  coefficient  at  the  source  and  receiving  ends.  These  relations  also  indicate  that  the 
difference  in  modal  velocities  of  propagation  can  lead  to  abrupt  changes  in  the  signal 
waveforms.  This  is  best  illustrated  in  Fig.  3.2  where  only  line  1  is  excited.  For  a  reading  at 
x  =0  on  line  2,  the  reflected  even  and  odd  mode  signal  (5  2  and  S'  2)  arrive  at  slightly 


different  times  and  since  their  contributions  are  of  opposite  signs  (cf.  Eq.(3.13a)),  they  pro¬ 
duce  an  impulse  in  the  waveform  of  line  2.  Such  a  disturbance,  however,  does  net  occur  in 
line  1  since  S  2  and  S'  2  are  added  according  to  (3.12a). 

Experimental  readings  were  performed  using  a  PG502  Tektronix  pulse  generator  and  a 
microstrip  pair  (Figs.  3.2  and  3.3).  The  readings  were  found  to  verify  the  theory  which  can 
also  be  extended  to  describe  complex  terminations. 

3.2  Capacitive  Terminations 

Most  digital  applications  for  transmission  lines  include  devices  such  as  diodes,  transistors, 
and  logic  gates.  Ideally  the  input  impedance  of  such  devices  is  very  high  and  the  input 
capacitance  is  minimized  to  insure  small  switching  delays.  In  practice  as  rise  and  fall  times 
become  shorter,  any  small  capacitance  must  be  taken  into  consideration,  this  is  more  impor¬ 
tant  when  several  identical  devices  are  connected  to  the  far  end  of  a  single  microstnp  line. 
Capacitance  is  related  to  the  carriers  in  the  devices  ;  for  instance,  in  the  case  of  a  bipolar 
transistor,  a  diffusion  capacitance  associated  with  the  charge  store  in  the  base  and  a  transi¬ 
tion  capacitance  related  to  the  base-emitter  space  charge  layer  make  up  the  total  input  capa¬ 
citance.  This  suggests  that  the  impedance  seen  is  strongly  nonlinear;  however,  an  average 
capacitance  can  always  be  defined  and  used  to  implement  a  good  linear  model. 

In  the  case  where  the  load  impedance  consists  of  a  resistor  R ,  in  parallel  with  a  capacitor 
C,  (Fig.  3.4a),  the  reflection  coefficients  T Le  and  F  u  for  the  even  and  odd  modes  in  the  fre¬ 
quency  domain  become  functions  of  frequency. 
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Figure  3.2  :  Cross-coupling  for  passive  terminations  with  far  end  opened.  Theoretical  (plots) 
and  actual  (photographs)  waveforms  for  a  long  (left)  and  short  (right)  pair  of 
coupled  microstrip  lines.  Photographs  :  top  waveforms  are  V ,(())  (driving 
line)  with  vertical  :  2  V  /div .  bottom  waveforms  are  V/0)  (idle  line)  with 
vertical  :  0.5  V  /div.  Each  horizontal  division  is  5  ns.'  The  spikes  on  the 
waveforms  of  V/0)  are  due  to  the  delay  differences  between  odd  and  even 
mode  reflections.  For  a  short  pair  of  lines,  these  spikes  can  no  longer  be 

detected  by  the  scope  because  of  their  short  duration. 


V2  (o)  (volts) 


Figure  3.3  :  Cross-coupling  for  passive  terminations  with  far  end  shorted.  Theoretical  (plot) 
and  actual  (photograph)  waveforms.  Photographs  :  top  waveform  is  V  ,(0) 
(driving  line)  with  vertical  ;  I  V  /div,  bottom  waveform  is  V/0)  (idle  or 
sense  line)  with  vertical  :  0.5  V  /div.  Bach  horizontal  division  is  5  ns . 
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TId  = 


~Pd  +  juTj 

1  +  jcOTj 


(3.31) 


where  we  have  defined 


R  -Zc 
R  +  Z„ 


RZ„C 
R  +  Z, 


(3.32) 


R  -Zd 
R  +Zd 


RZdC 
R  +  Zd 


(3.33) 


Since  T Le  and  T Ld  are  functions  of  frequency,  the  time  domain  solutions  for  the  modal 
voltages  are  convolutions  of  the  delayed  source  voltages  and  the  inverse  transforms  of 
Tu(cj)  or  T id(ai). 


ari-vMv-‘u 


-  £_+  2kl.)  +  Vj2(r  - 


(3.34a) 


Sfrjc)  -  —■  Z,.-y\  Iltv,  ,(t  +  - — 20+11)  +  VJt  +  i — * 

2\Z.e  +  i  — Q  ve  Vc 


F~  [rf  ♦*(«)] 


(3.34b) 


S'U,X)  2(Zd  +  Zs )  <§Vs  l(‘ 


_  *  +2kl)-  v  jt  _  £  +-2ft)]r^* 


F  [r Uu)] 


(3.35a) 


S’X:j)  =  -~  ■  £[Vsl(r  +  - — v  +  £ - M±1 )l  )r»  * 

2 \Zd  +  Zs)  tsQ  Vj  vd 


F  [n/'M 


(3.35b) 


where  *  denotes  a  convolution.  Let  us  consider 
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t  -1  1  +  n  y  +  1  ~ 

8(t )  +  (-1)* +l  - -  <?  f  (3.40b) 

j- 0  Tf 


^-7-Zl  7  \  Elv./I  -^2iL)-vj,  -  ±±2iL)iri,* 


*-i  l  +  \J+i  ~ 

SlJ  +  C-l^fly;  t'  — —  e  - 

/=0  Td 


(3.41a) 


y/l  pc )  =  — Elv.,0  +- — ?(*-+ — )  -Vt£t  +  - — 2(k  t  ?)z-]r/rf* 

**&d  +  A  >  k  =0  vd  vd 

*-l  1  +  n  ■/  +  1  71 

S(r)  +  (-iy+1I]D,i<+1C''  - -  e  *  (3.41b) 

j= O  Trf 

In  the  case  where  the  applied  excitation  is  a  pulse  (cf.  Fig.  3.4b),  V,  ,(r )  and  Vs  Jt )  have  the 


Vt(f)  = 


,  (j  ->.) 


r  (r  ~f  J 

p(h-0  ’ 


t  ,<{  , 


t2^t 


t  }^t 


(3.42) 


Moreover,  the  convolution  integral 


I- I 


( t  -  1 1— A)  -7 

T7^TTV- e  x'i* 


(3.43) 


can  be  solved  by  using  the  j’h  order  derivative  Theorem  for  Laplace  transforms. 
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A'  e~sKd A  = 


-fir  -  £-?rr  e' 


c  ./'  +  *  c  m  + 1 

°  m  =0  ■> 


Therefore,  if  we  define 


C=vsitf 


*  -j 

L±£. 

y+i 

8(f)  +  (-1>  'ED.*  t> 

j=0 

T 

e  T 


we  have 
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1  +p 

y  +  ! 

v 

T 
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k  -1 


C  =  Vs(r )  +  (-!)* '‘Hz?,,  r 


y=0 


1  +  P 


|y  +  i 


(3.44) 


(3.45) 


(3.46) 
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-(-1  y-'ZD;,  t' 

j= 0 


(t-t  j) 

P(t  3-0 
)/+i 

1+P 


j!  ry  +  1  -  £m!  rm+l  e 


t  r  — r  3 » 


m  =0 


vp 

3~l  4) 


1  f  -r  3) 

(j+l)!r>+2-  £m.'rm+1  (/-r^"m+1e  T 

m  -0 


for  f  3^f  <t4 


where  we  have  made  use  of 


/(f)*  80)0"*'  =  /(f)  (3.47) 

Using  these  relations  in  (3.40)  and  (3.41),  we  can  obtain  the  expressions  for  the  modal 
voltages.  These  equations  can  be  easily  implemented  on  a  computer  program  to  simulate  the 
waveforms  produced  by  a  pulse.  Such  example  is  shown  in  Figs.  3.5  and  3.6.  As  expected  the 
pulse  width,  the  length  of  the  lines,  and  the  RC  time  constant  are  the  most  important 
parameters  that  determine  the  shape  of  the  waveforms.  As  for  the  case  of  passive  loads,  fast 
pulses  can  result  at  the  sending  end  of  a  quiet  line,  (waveshape  of  V  /O)  ).  These  fluctua¬ 
tions  cannot  be  observed  accurately  using  presently  available  scopes  since  their  duration  is 
much  shorter  than  the  response  of  th|se  instruments. 


V2(o)  (volts) 


Figure  3.6  :  Cross-coupling  for  capacitive  loading.  Theoretical  (plot)  and  actual  (photograph) 
waveform  of  V  40)  (idle  line).  Photograph  :  voltage  reading  is  divided  bv  10. 
R  is  infinite  (open). 
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CHAPTER  4 


ANALYSIS  OF  MULTIPLE  LINE  STRUCTURES 

The  previous  chapters  dealt  with  with  two  microstrip  line  structures  for  which  methods 
of  measurements  and  time  domain  behavior  were  investigated.  Structures  consisting  of 
more  than  two  microstrip  lines  are,  however,  of  greater  occurrence  in  practical  situations.  In 
particular,  the  behavior  of  three-line  systems  in  the  frequency  domain  has  been  explored  by 
several  authors.  Collier  and  El  Deeb  [29]  have  determined  the  scattering  parameters  for  a 
six-port  refectometer;  Tripathi  [30]  derived  expressions  for  the  immitance  parameters  of 
symmetrical  three-line  microstrip  circuits.  Pavlidis  and  Hartnagel  [31]  derived  the  funda¬ 
mental  modes  of  propagations  for  these  structures.  Paul  [32]  constructed  the  matrix  chain 
parameters  for  multiconductor  transmission  lines  in  the  frequency  domain.  Pulse  propaga¬ 
tion  however,  requires  a  different  and  more  complex  analysis  in  the  time  domain.  For¬ 
tunately,  by  using  the  insight  gained  in  solving  the  two-line  problems,  expressions  for  mul¬ 
tiple  line  structures  in  the  time  domain  can  be  written  by  inspection  provided  that  the 
different  modes  of  propagation  are  identified.  The  goal  of  this  chapter  is  to  illustrate  such  an 
approach  and  generalize  the  techniques  for  n-line  systems. 

4.1  Three-Line  Structures  and  Modes  of  Propagation 

We  first  assume  that  the  three  microstrip  lines  shown  in  Fig.  4.1  are  identical  and  that 
the  two  edge  spacings  are  equivalent  ;  then  the  differential  equations  relating  line  voltages 
and  currents  can  be  written  directly. 


dVj  _  .  3/  , 

“  ^  II 

dx  Qt 


+  Ll2 


9/ 2 
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(4.1a) 
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_  av,  a/.  ,  a/ 2  .  .  a/3 

^  12  +  i-11 -  +  i'  12 - 

bx  y  y  y 


_E±~  ,  a;  i  .  ,  a/  2  ^  ,  a;  3 

~  L  13  +  i'  12 -  +  L  || - 

a*  y  y  y 


(4.1b) 


(4.1c) 


_ iii _r  evi  d^2  av3 

—  <-11  +1-12 -  +  C  |i- - 

a*  at  at  a? 


(4.2a) 


§L±-r  £Ii+r  ^  j.  r  9V> 
<-  12 -  +  <-  || -  +  C  - 

a*  a?  y  y 


(4.2b) 


§L±-r  av,  av,  av3 

-  O  13  +<-12 -  +C|1 - 

a*  y  y  at 


(4.2c) 


where  the  LtJ  and  the  Ct]  are  related  to  the  physical  parameters  by 


L  a  =  L, 


L  12  —  Lm 


L  i3  —  Ln 


c  „  =  cs  +  cm  +  c„ 


c  12  =  — C„ 


c  i3  =  C„ 


The  subscript  m  denotes  mutual  parameters  between  adjacent  lines  and  the  subscript  n. 
mutual  parameters  between  nonadjacent  lines.  If  we  define 
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Va  =  V,-V3 


(4.9a) 


1=J>-1 


(4.9b) 


we  can  determine  a  first  mode  of  the  system  by  subtracting  (4.1c)  from  (4.1a)  and  (4.2c) 
from  (4.2a).  This  yields 


—  -Un-W¬ 

ax  y 


(4.10a) 


d*  Qt 


(4.10b) 


The  velocity  of  propagation  associated  with  this  mode  is 


/(L  U~L  l3*C  U-C  ti)  V(L,  XCS  +C,„  +2CJ 


(4.11) 


Its  characteristic  impedance  is 


'L  i)  —  L  i3 

('ll  ~  C  13 


'  ^~Ln 

C,  +  Cm  +  2C„ 


(4.12) 


A  second  mode  can  be  obtained  by  defining 


V(  =  V,+fV2+  V: 


(4.13a) 


/,-  =  /,  +  V  2  +  /  ; 


(4.13b) 


where  £  is  a  constant  to  be  determined.  When  this  linear  combination  is  made  using  (4.1) 
and  (4.2),  we  get 


-  + fL ,2  + L,,)— 

6*  9 1 


*  s  a* 


+  U*13  +  O- .2  +  0-1.)— 

and 


=  (CU  +  |CI2  +  C13)^- 


+  (^  +  c„  +  %M^ 

*  £  # 

+  (C,3  +  |C12  +  C„)^ 

8* 

We  then  introduce  the  approximations 


£  u  +  O-  12  +  L  i3  **  £  ii  +  O  l2 

^  n  +  i2  +  c  13  ^  c  ii  +  £c  i2 

which  reduce  Equation  (4.14)  to 


~  (L“  +  +  (Lll+2-^)£'~  +  (Lu  +  ^l2)- 


e/ 


Ox  =  (cn  +  ^Cl2)‘^i  +  (cm  +  2'^li)£~i  +  (cn  +  £c,2) 


If  we  choose  £=+>/2,  we  then  have 
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(4.18) 


5— =(L„  +  x/iLp)-^-1 
dX  Qt 


2ii-=(Cu  +  ^Cp)^- 


(4.19) 


This  pair  of  equations  defines  a  second  mode  with  propagation  velocity 


/(L  ,,  +  x/2L  pXC  ,i  +  x/2C12)  Va,  +  'f2LmXCs  +  C„  +  (l-V2)Cm) 


(4.20) 


and  an  associated  wave  impedance 


(I,i  +  412) 


V  (Cu  +  n/^C,,) 


(Ls  4-  V2Lm) 

(Cs  +  C„  +(l-x/2)Cm) 


(4.21) 


Finally,  a  third  mode  is  obtained  in  a  similar  manner  by  defining 


V„=  V,+T)V,+  V; 


(4.22a) 


I  i)  =  I  I  +  Tjl  2  +  I . 


(4.22b) 


Making  the  same  approximations  as  per  Equations  (4.15)  and  (4.16)  and  setting  ?)=—'/ 2,  we 


^L=(Lu_^)^L 

6*  3? 


(4.23a) 


— =  (C —  x/2)—- — 


(4.23b) 


The  propagation  velocity  is 
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Va  „  -  y/2 L  ,,XC „  -  'he ,,)  V(Ij  -  v/2Lm  XCS  +  Cm  +(l+V2)Cm) 


(4.24) 


and  the  characteristic  impedance  is 


rnnL 

V  c„-V2cl2  V  cs  +  c„  +  (i+v 


(l+>/2)C„ 


The  line  variables  can  then  be  expressed  in  terms  of  the  three  modal  variables  as 


(4.25) 


V,  =  - 


v£  +  v„ 

V  +  - 1 


(4.26) 


(4.27) 


V,=  4- 


LV 


(4.28) 


and 


/i 

/2 


I  +  n 


Ij-In 

2V2 


'3 


I  n  +  I  n 
2 


(4.29) 

(4.30) 

(4.31) 


It  must  be  emphasized  that  these  results  apply  if  Equations  (4.15)  and  (4.16)  are  valid.  In 
general,  such  approximations  are  very  well  justified  since  in  most  cases 
C„  «C, +(l  ±  V2)Cm  and  L„  «LS  ±'JlLm.  Figure  4.2  shows  the  different  excitations 
associated  with  the  three  modes  of  the  system  in  the  case  where  the  mutual  coupling 
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where 
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\  o 


Z«zl± 

Za  +  2, 


z.-zz 

za  +  ZL 


(4.35) 


r 


si  ~ 


z£-z, 
zi  +  zs 


ru  = 


z£-z, 
ZL  +  zL 


(4.36) 


_  z,-z5  _  zv  —  Z i 

71  Z„  +  Z5  iT>  Z^  +  Zi 


(4.37) 


The  modal  voltages  are  as  defined  by  Equations  (4.9),  (4.13),  and  (4.22)  where  the  subscript 
/  denotes  a  forward  moving  wave  and  b,  a  backward  moving  wave.  Using  (4.26M4.31), 
the  line  variables  can  be  obtained.  If  we  apply  the  boundary  conditions  at  x  =0  and  at  x  -l , 
we  get 


(4.38) 

(4.39) 

(4.40) 


4.2  N-Line  Structures  and  Generalization 

From  the  analysis  performed,  we  can  anticipate  that  n-line  systems  have  solutions 
analogous  to  three-line  structures.  More  precisely,  if  the  terminations  are  identical,  an 
expression  can  be  written  for  the  n-line  structure  from  the  modal  characteristics.  Since  a 
matrix  representation  is  more  convenient,  the  problem  involves  the  determination  of  the 
eigenvalues  and  eigenvectors  which  can  be  obtained  using  some  approximations.  The  Tele¬ 
graph  Equations  then  become 


where  L  and  C  are  nxn  matrices  and  V  a  nxl  matrix  for  the  line  voltages.  If  the  n 
modes  of  the  system  can  be  found,  we  can  define  an  eigenvector  matrix  E  such  that 


V  m  =  E  V  (4.42) 

V m  is  the  column  matrix  for  the  n  modal  voltages.  The  solution  in  terms  of  E  and  the 
matrix  associated  with  the  voltage  sources  at  x  =0,  V  f,  is 


v fm  =  Am  YjE  rsmrlm  v,(ufmt  (4.43a) 

k  =0 

v („„  =  -Am  ££  r,mrLm  v bm(ubm)  (4.43b) 

k  =0 

where  V  frn  and  V bm  are  the  modal  voltage  matrices  for  the  forward  and  backward  waves, 
respectively.  A,„,  Tj,,,,  and  are  diagonal  n  xn  matrices  associated  with  the  source  and 
load  impedances  and  the  corresponding  mode  ;  ufm  and  uhm  are  the  arguments  associated 
with  the  mode  m  for  forward  and  backward  variables,  respectively 


'■fm 


x  +  2  kl 


Abrn 


=  t  + 


x  -  2(*+i); 


(4.44) 


where  vm  is  the  propagation  velocity  for  mode  m.  From  these  relations,  the  line  voltages 
can  be  found  by  applying 


V  =£ 


Vfm  +V, 


(4.45) 


We  then  observe  that  the  problem  of  finding  the  signal  magnitudes  at  any  time  or  any 
position  of  an  n  -line  structure  becomes  that  of  finding  the  eigenvalues  and  the  eigenvector 
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associated  with  the  system.  For  a  large  number  of  lines,  this  problem  becomes  nontrivial 
however,  by  using  computer  routines  or  by  setting  some  approximations  such  as  those  in 
Equations  (4.15M4.16),  one  can  arrive  at  satisfactory  results.  For  instance,  it  can  be 
assumed  that  mutual  parameters  associated  with  any  two  nonadjacent  lines  are  negligible. 
This  generates  L  and  C  matrices  with  only  principal  and  secondary  diagonal  elements  and 
facilitates  the  analytical  or  computer  task. 

The  solution  to  n-line  pulse  driven  structures  applies  to  many  digital  network  problems. 
When  n  logic  gates  are  switching  n  other  gates  via  transmission  lines,  noise  and  reflections 
induced  on  a  quiet  line  can  be  determined  at  any  point.  This  would  provide  margins  for  the 
applied  signals  (magnitude,  pulse  width,  rise  and  fall  times),  and  the  microstrip  line 
geometries  (spacing,  width,  and  dielectric  constant). 


CHAPTER  5 


CONCLUSION 

This  study  explored  some  of  the  problems  associated  with  two-line  structures  and 
attempted  a  generalization  to  multiple  line  systems.  These  analyses  represent  the  first  steps 
for  the  investigation  of  more  complex  situations  involving  nonlinearities,  capacitive  termi¬ 
nations,  and  discontinuities  along  the  lines.  N-line  systems  with  arbitrary  terminations  can 
therefore  be  treated  using  a  numerical  approach  to  provide  reliable  design  guidelines  for 
digital  networks. 
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